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ABSTRACT 

In this paper, Jupiter's Great Red Spot (GRS) is used to determine properties of the Jovian 
atmosphere that cannot otherwise be found. These properties include the potential vorticity of 
the GRS and its neighboring jet streams, the shear imposed on the GRS by the jet streams, and 
the vertical entropy gradient (i.e., Rossby deformation radius). The cloud cover of the GRS, 
which is often used to define the GRS's area and aspect ratio, is found to differ significantly 
from the region of the GRS's potential vorticity anomaly. The westward-going jet stream to 
the north of the GRS and the eastward-going jet stream to its south are each found to have 
a large potential vorticity "jump". The jumps have opposite sign and as a consequence of 
their interaction with the GRS, the shear imposed on the GRS is reduced. The east-west 
to north-south aspect ratio of the GRS's potential vorticity anomaly depends on the ratio of 
the imposed shear to the strength of the anomaly. The aspect ratio is found to be ?a2:l, but 
without the opposing jumps it would be much greater The GRS's high-speed collar and 
quiescent interior require that the potential vorticity in the interior be approximately half that 
in the collar. No other persistent geophysical vortex has a significant minimum of potential 
vorticity in its interior and laboratory vortices with such a minimum are unstable. 



1. Introduction 

llonly a relatively thin (~10 km) outer layer of Jupiter's atmosphere containing the 
visible clouds and vortices is accessible by direct observation. Most of the details of 
the underlying layers, such as the vertical stratification, must therefore be determined 
indirectly. In this paper, we present one such indirect method. In particular, we use 
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the observed velocity field of a persistent Jovian vortex to determine quantities rele- 
vant to both outer and underlying layers. These quantities are the potential vorticity 
of the vortex, the potential vorticity of the neighboring jet streams, the flow in the 
underlying layers, and the Rossby deformation radius L^, which is a measure of the 
vertical stratification. We demonstrate the method using Voyager 1 observations of 
the Great Red Spot (GRS). 

We are not the first to use the GRS velocity field as a p robe of the Jovian atmo- 
sphere (IDowhng and Ingersolllll988l . ll989l : ICho et al.ll200lb . However, our approach 
differs from previous ones in several significant respects. First, the GRS velocity 
field is sufficiently noisy that we do not, unlike in previous analyses, take spatial 
derivatives of the velocity to compute potential vorticity. Instead, we solve the in- 
verse problem: We identify several "traits" of the GRS velocity field, where a trait 
is a feature of the velocity field that is unambiguously quantifiable from the noisy 
data. We then construct a model for the flow and determine "best-fit" values for the 
model parameters such that the model velocity field reproduces the observed traits. 
Furthermore, for a given set of parameter values, we construct the model velocity 
field so that it is an exact steady solution of the equations that govern the flow. For 
the Voyager 1 data, we find that a best-fit model (i.e., a trait-reproducing steady so- 
lution) determined in this manner agrees with the entire GRS velocity field to within 
the observational uncertainties. 

A second way in which our study differs from previous ones is that we explicitly 
compute the interaction between the GRS and its neighboring jet streams. We show 
that the interaction controls the aspect ratio of the GRS's potential vorticity anomaly, 
which is relevant to recent obs ervations that show t he aspect ratio of the GRS's cloud 
cover to be a function of time (^Simo n-Miller et aL 2002,) . The changing cloud cover, 
if symptomatic of changes in the GRS's potential vorticity anomaly, would be indica- 
tive of a change in the interaction and a corresponding change in the best-fit values of 
the parameters that govern the interaction. Finally, in this study, we quantify the re- 
lationship between individual traits and individual parameters. When a trait is nearly 
independent of all parameters except for one or two, a clear physical understanding 
is obtained between "cause" (a model parameter) and "effect" (a GRS trait). 

Our philosophy is to use a model with the fewest free parameters that is an exact 
steady solution to the least complex governing equation, yet can still reproduce the 
observed velocity to within its uncertainties. The danger of more complex models 
is that they have larger degrees of freedom. By varying parameters they can fit the 
observed velocity but misidentify the relevant physics. For the Voyager 1 data con- 
sidered here, we use the 1.5-layer reduced gravity quasigeostrophic (QG) equations 
and a model with nine free parameters. The Voyager 1 data can be reproduced with 
this model and does not warrant models with more free parameters or governing 
equations with more complexity. 

The rest of the paper is organized as follows. In §2 we determine the GRS veloc- 
ity field from Voyager 1 observations and then identify traits of the velocity. In §3 we 
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review the governing equations and describe a decomposition of the flow around the 
GRS into a near-field, a far-field, and an interaction-field. In §4 we define the model 
and list its free parameters. In §5 we determine best-fit parameter values, i.e., param- 
eter values for which the model reproduces the traits. In §6 we discuss the physical 
implications of the best-fit model, and in §7 conclude with an outline for future work. 

2. GRS velocity field 



a. Determination of GRS velocity 

In [Mitchell et al.l (Il98lh . Voyager 1 images were used to determine the GRS ve- 
locity field by dividing the displacement of a cloud feature in a pair of images by 
the time interval between the images (typically one Jovian day or sslO hours). The 
cloud features were identified by hand rather than by an automated approach such as 
Correlation-Image- Velocimetry (CIV: Fincham and Spedding 1997), and may there- 
fore contain spurious velocities on account of misidentifications. Furthermore, di- 
viding cloud displacement by time does not account for the curvature of a cloud 
trajectory, since in ten hours, a cloud feature in the high-speed collar travels almost 
a third of the way across the GRS. However, due to the unavailability of the original 
navigated images, we use the Mitchell velocities, but remove some of the errors by 
a procedure described in appendix A. The procedure leads to the removal of 220 of 
the original 1100 measured cloud displacements and the addition of 7100 synthetic 
measurements. The net result is that the uncertainty in the velocity field is reduced 
from « 9 m s^^ to « 7 m s^^ Fig.HJshows the processed GRS velocity field. Consis- 
tent with previous analyses, the velocity field shows a quiescent core and high-speed 
collar. The inner part of the collar has anticyclonic vorticity, and the outer part has 
cyclonic vorticity. The peak velocities in the collar are 140 it 7 m s^^ and the peak 
velocities in the core are 7 it 7 m s^^ 

The GRS is embedded in a zonal (east-west) flow. The zonal mean of this flow, 



avera ged over 142 Jovian days, was computed from Voyager 2 images (ILimaye 



1986 ). and is shown in Fig. [2l Between 15°sQ and 30°S, the profile is character- 
ized by a westward-going jet stream that peaks at w 19.5°S, and an eastward-going 
jet stream that peaks at ?a 26.5°S. The uncertainty in t he profile is 7 m s^^ (Note that 



most likely due to navigational errors (|Limaydl 19861 ). the published profile must be 



shifted north by 0.5° so as to be consistent with the navigated latitudes of Voyager 1 
in Fig. [U) The GRS was observed to drift w estward at a rate of 3-4 m s ~^ with 



respect to System III during the Voyager epoch (|Dowling and Ingersolllll988h . 



b. Pitfalls to be avoided when analyzing GRS velocity 

We do not compute quantities by taking spatial derivatives of the velocity data, as 



In this paper all latitudes are planetographic. 
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Fig. 1. Velocity vectors of the GRS with respect to System III as determined from Voyager 1 images. 
The velocities were determined by d ividing the displacem ent of a cloud feature in a pair of images 

' letall m 



by the time between the two images ([Mitchell 



19811) . and then correcting for the fact that cloud 



trajectories over typical image separation times are not straight lines (see §2a). The 23°S latitude is 
defined to be the principal east-west (E-W) axis, and the 77°W longitude is defined to be the principal 
north-south (N-S) axis. 



this tends to amplify small length scale noise. For example, in the high-speed collar, 
we found that the uncertainty in vorticity obtained by differentiating the velocity is 
K, 35% of the maximum vorticity. If vorticities must be found, it is usually better 
to integrate the velocity to obtain a circulation and then divide by an area to obtain 
a local average vorticity. We also do not average the velocity locally, which is a 
standard way of reducing noise. For example, if the GRS velocity is averaged over 
length scales greater than 2°, the peak velocities and vorticities are severely dimin- 
ished. This is due to the fact that an averaging length of 2° is too large; it corresponds 
to K, 2500 km, which is the length scale over which the velocity changes by order 
unity (cf. the width of the high-speed collar). Finally, we do not obtain a quantity 
by adding two numbers of similar magnitude but opposite sign, so that the resulting 
sum is of order or smaller than the uncertainty in each of the numbers being summed. 
For example, if the velocity is assumed to be divergence-free, the vertical derivative 
of the vertical velocity 3vj/3z can be obtained by computing the negative of the hor- 
izontal divergence 3v_r/3.x + 3vy/33' . However, a simple scaling argument shows that 
the horizontal divergence is smaller than each partial derivative term separately, and 
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Fig. 2. Averaged zonal velocity v°° in System III jLimavdl 19861) 



in particular, is of the same order as the uncertainty in each term (which is relatively 
large because the terms are derivatives of noisy data). Thus dv,/dz computed in this 
fashion would have order unity uncertainties. 



c. Traits of GRS velocity 



The traits that we consider are derived from the north-south (N-S) velocity along 
the principal east-west (E-W) axis and from the east-west (E-W) velocity along the 
principal north-south (N-S) axis. The E-W and N-S principal axes are defined to be 
the 23° S latitude and the 77° W longitude respectively. The point of intersection of 
the principal axes is roughly the centroid of the GRS as inferred from its clouds. The 
velocity profiles along the axes are shown in Fig. [3] To better understand the pitfalls 
of local averaging, Figs. [3K and 3c show the velocities from Fig. 1 for points that lie 
within ±0.7° of the axes, while Figs. 3b and 3d show the velocities that lie within 
±1.4° of the axes. The axes labels x and y in the figure denote local E-W and N-S 
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cartesian coordinates. Based on the figure, we define tlie following to be traits of tfie 
velocity field: (1) the northward-going jet and southward-going jet in Figs. [3ta)-(b) 
that peak a.tx = ±9750 ± 500 km respectively and have peak mag nitude V^f^ = 95 ± 
7 m s^\ (2) the small magnitude N-S velocity in < 6000 km, (3) the eastward- 
going jet and westward-going jet in Figs. [3];c)-(d) that peak aty = —3500 ± 500 km 
and y = 5500 ± 500 km respectively, and have peak mag nitude 140±7 m s^^, (4) 
the small magnitude E-W velocity in |j| < 2000 km. Traits (2) and (4) illustrate the 
quiescent interior of the GRS. Traits (1) and (3) illustrate the high-speed collar. The 
uncertainties in peak velocities are from the global estimate in §2a. The uncertainties 
in peak locations are not rigorous. They are from an estimate of the spatial scatter of 
points near the peak location. Henceforth, traits (1) and (2) will be referred to as the 
N-S velocity traits, and traits (3) and (4) as the E-W velocity traits. 

3. Governing equations 



a. 1.5-layer reduced gravity QG approximation 

We do not model the whole sphere, but only a domain that extends from 15°S 
to 30° S. For the flow in t his domain, we adopt the 1.5-layer reduced gravity QG 



equations on a beta-plane d lngersoU and Cuong|1981h. A derivat ion of the equations 



and the justification for their use can be found in lDowlingI (119951) . Briefly, the layers 
correspond to an upper layer (also called "weather" layer) of constant density p i and 
a much deeper lower layer of constant density P2 > Pi- The upper layer contains the 
visible clouds and vortices while the lower layer contains a steady zonal flow. The 
two layers are dynamically equivalent to a single layer with rigid bottom topography 
hh and effective gravity g = g/(p2 — Pi)/P2, where gj is the true gravity in the weather 
layer, and the bottom topography is a parametrization of the flow in the lower layer. 
The governing equation for the system advectively conserves a potential vorticity q: 

^.(|.v.v),.0. 

^vv-;!^+§^?i+p,. (2) 

Here x and y are the local E-W and N-S coordinates, \|/ is the streamfunction, v = 
z X V\|/ is the weather layer velocity, z is the local vertical unit vector, P is the local 
gradient of the Coriolis parameter f{y), fo is the local value of f{y), and is the 
local Rossby deformation radius. Since g appears only in combination with ht,, we 
shall refer to ghb as the bottom topography. Restricting ght, to be a function of y alone 
restricts the flow in the lower layer to be steady and zonal with no vortices. The case 
ghh = 0, or "flat" bottom topography, corresponds to the lower layer being at rest in 
the rotating frame. 
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Fig. 3. Parts (a) and (b) show the north-south component of the velocity from Fig. [T] for points that 
lie within 0.7° and 1.4° respectively, of the principal east-west axis. Parts (c) and (d) show the east- 
west component of the velocity for points that lie within 0.7° and 1.4° respectively, of the principal 
north-south axis. Parts (a), (b), (c), and (d) contain 903, 1992, 586, and 1163 points respectively. 



b. The near-field 

We assume that the GRS is a compact region (or patch) of anomalous potential 
vorticity. We denote the potential vorticity distribution of the GRS by q^^^{x,y), and 
for reasons that will become clear below, we refer to q^^^ as the near-field. We define 
the streamf unction and velocity of the near-field to be: 

^g«^(x,3;)^(v2-1/l2)x)/^«^(x,3^) (3) 

v'^^'^ = z X WMf^^^. (4) 



The velocity induced by a QG patch decays as exp(— r/L^), where r is the distance 
from the patch boundary (IMarcusI Il990l Il993l) . Due to the exponential decay of 
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velocity, a region of fluid that contains the patch and whose average radius is a few 
Lr greater than the patch radius will have a circulation (or integrated vorticity) that 
is approximately zero. It would therefore be incorrect under the QG approximation 
to refer to the vorticity of the GRS as anticyclonic since its net vorticity is zero. On 
the other hand, the potential vorticity of the GRS is anticyclonic, as is the vorticity 
of most of its quiescent interior and the inner portion of its high-speed collar, but 
the vorticity of the outer portion of its collar is cyclonic (which is easily verified by 
noting that the azimuthal velocity in that region falls off faster than 1/r). 

c. The far-field 

The region of flow two or three deformation radii distant from the patch boundary, 
where the influence of the GRS is small, is defined to be the far-field. We assume the 
far-field flow to be zonal and independent of time and longitude. Eq.[2]then provides 
a relationship between the far-field velocity v°° = v~(y)x, and the far-field potential 
vorticity q°°{y): 

For all calculations in this paper, v'^{y) is prescribed from Fig. |2] and the correspond- 
ing streamfunction \|/°° from — / v" dy. At 22.5°S, which is the center of the domain, 
P = 4.6 X 10"^2 and /o = -1.4 x 10""^ s^^ Thus if were known, Eq.[5] 

shows that specifying q'" is equivalent to specifying ghi,. 

d. The interaction field 

Let q{x,y) be a steady solution of QG Eqs. [T]-[2] that consists of an anomalous 
patch of potential vorticity embedded in a far-field flow that is zonal and steady. We 
decompose q into three components: 

q{x,y) ^ q-{y)+q''''\x,y)+q"'^{x,y). (6) 

The superposition of q'^^^ and q°° is not an exact solution because the far-field flow 
is deflected around the patch. We define the interaction potential vorticity to 
represent the deflection of flow such that the total q given by Eq.[6]is an exact solution 
of the QG equations. Note that by definition, q^^^ asymptotes to zero both far from 
and near the patch. We define the interaction streamfunction and velocity to be: 

q^^\x,y)^{y'-llLl)^^'^\x,y) (7) 

v^^^ = z X VV^^. (8) 

With these definitions, the total velocity v and its streamfunction \|/ are superpositions 
of the near, interaction, and far-field components: 

\|/ = x|/°° + + \|/^^^ (9) 
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v = v~ + v^^^ + v^^'^. (10) 

Note that in the hnear relationships between the potential vorticity and streamfunc- 
tion given by Eqs. |5] [3l and |7J it is only the far-field component that contains the 
inhomogeneous bottom topography and p terms. 

4. Model definition 



a. Model for far-field q°° 

Laboratory e xperiments j somm eria et aL |[l989|; ISolomon et al. I 1993b and numer- 
ical simulations ( Cho and Po lvani 1996l: lMarcus et al" 20od) show that if the weather 
layer is stirred and allowed to come to equilibrium, the potential vorticity organizes 
itself into a system of east- west bands. The bands have approximately uniform q and 
are separated by steep meridional gradients of q. The meridional gradients are all 
positive (i.e., have the sign as P) so that q{y) monotonically increases from the south 
to the north pole like a "staircase". The corresponding v°° has alternating eastward- 
going and westward-going jet streams, with eastw ard-going jet stream s occurring at 
every gradient or "jump". Recent measurements (|Read et al.l l2006bh of the Jovian 
q^ are not entirely consistent with this picture, for they show gradients near both 
eastward-going and westward-going jet streams. We therefore model <7°° between 
30°S and 15°S by: 



i=\ 



tanh 



5,- 



+ 1 



(11) 



The jumps for this model occur at yi and have strength AQ, , where AQ,- can be positive 
or negative. The strictly positive 5; are a measure of the steepness of each jump. For 
all results in this paper, the jump locations were fixed at = 26.0°S and y2 = 20.0°S, 
which are near the jet streams in Fig. [2l The free parameters for the model are AQ, 
and 5; for i=l,2. (Models with up-to four jumps near each jet stream were also tested, 
and the results were consistent with the ones presented here.) 



b. Model for near-field q^'^^ 

We model the spatially compact q^^^ as a piecewise-constant function obtained 
by the superposition of M nested patches of uniform potential vorticity. The patches 
are labelled / = 1,2,- •• ,M from innermost to outermost patch. The principal E- 
W diameter of a patch is denoted by (Dv),-, and we define qf^^ such that q^'^^ = 
q^^^ within the boundary of the innermost patch (/ = 1), ^^^^ = q'^^^ between the 
boundary of the innermost patch (/ = 1) and the boundary of the next larger patch 
(/ = 2), and so on. The free parameters for the model are M, qf^^, and {Dx)i for 
= 1,2, ••• ,M. Once the free parameters for q'^'^^ and <7°° are specified, along with 
the value of L^, the iterative method given in appendix B can be used to compute the 
interaction-field such that the total ^ is a steady solution of the governing equations. 
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Table 1 . Relationship between an observable trait of the Great Red Spot and the model parameter that it 
determines. The ordering of the table, from top to bottom, is the order in which each model parameter 
is determined. A rigorous justification for the ordering is given in appendix C. 



Observable trait 


Model Parameter 


Distance between peaks in N-S velocity along E-W axis 
Magnitude of peak N-S velocity along E-W axis 
N-S velocity along E-W axis for x > {Dx)\/2 
N-S velocity along E-W axis for x < {Dx)\/2 
E-W velocity along N-S axis 


E-W diameter of GRS's potential vorticity: {Dx)\ 

Family of possible q^^^ and L,- 

Unique q^^^ and from family 

GRS's interior potential vorticity: q2^^, {Dx)2 

Far-field potential vorticity: AQi, 6, 



Note that the shapes of patch boundaries are not free but are also computed by the 
iterative method. 

5. Determination of best-fit parameter values 

a. Decoupling of N-S velocity traits from far-field q°° 

Here we show that the N-S velocity traits are insensitive to the far-field potential 
vorticity described by Eq. [TT] Fig. |4] shows a model computed using the iterative 
method in appendix B for the parameter values given in Table |2l The middle col- 
umn of Fig. |4] shows that for the E-W velocity along the N-S axis, all three velocity 
components, v°°, v*^^"^, and v^^^, contribute significantly. However, the rightmost 
column shows that for the N-S velocity along the E-W axis, \°° has no contribution 
(by definition), and the contribution of v^'^^ is negligibly smalH. Only v^^^ con- 
tributes significantly. We therefore conclude that the N-S traits depend primarily on 
parameters associated with q^'^^ and are insensitive to parameters associated with q°°. 
This decoupling leads to a logical order for determining the best-fit parameter values. 
The ordering is given in Table [T] and begins with the determination of (Av)i from the 
N-S traits. A more rigorous justification for the ordering is given in appendix C. 

b. Determination of best-fit Lr and q^^^ from N-S velocity traits 

Here we show that an M = 2 model is sufficient to capture the N-S velocity traits 
to within the observational uncertainties. For brevity, the terms interior and exterior 
are used in reference to the regions |x| < (Dv)i/2 and |x| > {Dx)i/2 respectively. 

For M = 1 , q^^^ is a patch of uniform potential vorticity. Models were computed 

^ The contribution of v"*'"^ is small because q'^^ comprises two highly (E-W)-elongated slivers 
north and south of the GRS (first column, bottom row of Fig.|4}. The associated v^'*'"^ follows highly 
(E-W)-elongated closed streamlines approximately concentric to q'^^ . Therefore, along the E-W axis, 
v^'^^ is primarily in the E-W direction. 
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Fig. 4. (left to right) Potential vorticity, E-W velocity along the N-S axis, and N-S velocity along the 
E-W axis. A gray-scale map is used for potential vorticity with black representing the most cyclonic 
fluid and white the most anticyclonic fluid, (top to bottom) The components due to the total q, the 
components due to alone, the components due to q'-'^^ alone, and the components due to q'^^ 
alone. The total ^ is a (uniformly translating) solution of the QG Equations. It was computed for the 
parameter values in Table|2]using the iterative method in appendix B. Note that the slivers of q^^^ (left 
column bottom row) are due to the displacement of the "jumps" or steep gradients in q°° as they follow 
streamlines that deflect around the GRS. 



for different values of (Dv)i, q'^^^ , and L^. For each model, the peaks of the N-S 
velocity along the E-W axis were found to occur ntx = ±{Dx)]_/2. The best-fit value 
of {Dx)\ = 19500 km was therefore inferred from trait 1 in §2c. Next, the best-fit 
values of 17^'^'^ and Lr were constrained using the observed peak magnitude V^^^ of 
the N-S velocity along the E-W axis. In particular, for a given value of L,-, the value 
of q^^^ was chosen so that the model reproduced the observed peak magnitude. By 
repeating this process for several values of Lr, a two-parameter family (i.e., a function 
of Lr and q^^^) of models that simultaneously capture the observed peak locations 
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Table 2. Best-fit parameter values for Lr, q'^^^, and q 



Parameter 


Best-fit value 


L, 


2400 km 


M 


2 




10.5 X 10"^ 


MRS 
H2 


6.0 X 10^5 


(O.v)l 


19500 km 




12000 km 




1.9 X 10^5 




-5.6 X 10"^ 


81 


300 km 


S2 


1000 km 



and peak magnitude was obtained. Some family members are shown in Fig [51 Note 
that the models do not capture the observed width of the northward and southward- 
going jets. In particular, for sufficiently small L,. {Lr = 1300 km), the model captures 
the rate of velocity fall-off in the interior but not in the exterior. For sufficiently large 
Lr {Lr = 2400 km) the opposite is true. For other values of L^, the rate of fall-off 
is too fast or too slow in both regions. To overcome this shortcoming, models with 
M = 2 were considered. 

For M = 2, q^'^^ is the superposition of two nested patches of uniform potential 
vorticity. Best-fit values of Lr, qf^^, and (Av)i were taken from the M = 1 model 
in Fig. [5] that captures the velocity profile in the exterior. As shown in Fig. [6l for 
^GRS _ qGRS^ jj^g M = 2 model does not capture the velocity profile in the interior. 
However, if ^^'^'^ is changed holding all other parameters fixed, only the interior flow 
changes (provided (Dy)i — (0 ^)2 > 2L^). Thus with all other parameters fixed, a 
genetic algorithm (IZohdill2003l) was used to determine values for q2^^ and (Ar)2 that 



minimize the L2-norm difference between the model N-S velocity along the E-W 
axis and the observed velocity in the interior. The parameter values obtained are 
listed in Tableland Fig.|7]shows that the N-S traits are captured for these parameter 
values. (Models with M = 3,4 were also tested, and the results were consistent with 
the ones presented here.) 

c. Determination of best- fit q"" from E-W velocity traits 



The best-fit q°° was determined from the E-W velocity traits. In particular, with 
the q*^^^ p arameters fixe d at their best-fit values from the preceding section, a genetic 
algorithm dZohdil 120031) was used to determine values for AQ, and 8; that minimize 
the L2-norm difference between the model E-W velocity along the N-S axis and the 
corresponding observed velocity. The parameter values obtained are listed in Table|2] 
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Fig. 5. N-S velocity along the E-W axis for models with M = 1. Crosses are Voyager 1 data from 
Fig.HJi. Solid black curve has qf^^ = 10.5 x 10"^ s"' and L,- = 2400 km. Dashed curve has qf^^ = 
6.5x 10"^ s"' andL,- = 3800km. Solid gray curve has qf'^^=l9.5 x 10"^ s"' and L,. = 1300 km. All 
models have the best-fit value of (Dx)i = 19500 km and were computed using the iterative method in 
appendix B. The N-S velocity along the E-W axis is insensitive to parameters of the far-field because 
of decoupling as described in §5a. 



Fig. [8] shows that the E-W traits are captured for these parameter values. The cor- 
responding q°° is shown in Fig. |9l The velocities for this trait-capturing model were 
found to match the GRS velocities in Fig. 1 to within the observational uncertainties. 

6. Physical implications of best-fit model 
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Fig. 6. N-S velocity along the E-W axis for models with M = 2. Each model has Lr, (£>a)i, qf^^, 
and (0^)2 set to the best-fit value in Table|2l but the values of ^2*'^ for each model are different. The 
dashed curve has ^2*'^ = g^*'^, the solid gray curve has ^2*'^ = O.S^j'^^ and the solid black curve has 
the best-fit value of ^j''^'^ = O.STq'^'^^. Crosses are Voyager 1 data from Fig. [3^. The N-S velocity along 
the E-W axis is insensitive to parameters of the far-field because of decoupling as described in §5a. 



a. Cloud morphology and GRS's potential vorticity anomaly 

The models show that the peak north-south velocities along the principal east- 
west axis occur at .x = ±{Dx)\/2, where {Dx)\ is the principal east- west diameter of 
the GRS's potential vorticity anomaly. Thus the best-fit value of {Dx)\ = 19500 km 
was inferred from trait 1 in §2c. In fact, the models show that not just the east- 
west extremeties, but the entire boundary of the GRS's potential vorticity anomaly 
is demarcated by the locations of peak velocity magnitude (|v|). This implies that 
an estimate of the area and aspect ratio of the GRS 's potential vorticity anomaly can 
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Fig. 7. Solid line is the N-S velocity along the E-W axis for a model with best-tit parameter values 
given in Table|2] Crosses are Voyager 1 data from Fig.[3^. 



be made directly from the observed velocity field without first determining a best-fit 
model. Traditionally, the clouds associated with the GRS have been used to infer 
the area and aspect ratio of the vortex. The east-west diameter of the cloud cover 
associated with the GRS is «24000-26000 km in length, which is «25% longer than 
the east-west extent of the potential vorticity anomaly as determined by our best-fit 
model. The north-south diameter of the cloud cover is also ss25% longer than that 
of the anomaly, so any estimate of the size of the GRS based on its cloud images 
rather than on its velocity overestimates the area of the potential vorticity anomaly 
by ^50%. 
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Fig. 8. Solid line is the E-W velocity along the N~S axis for a model with best-fit parameter values 
given in Table|2] Crosses are Voyager 1 data from Fig.[3^. 



b. Rossby deformation radius 



The models show the rate of fall-off of the north-south velocity in the outer por- 
tion of the high-speed collar to be almost independent of all parameters with the 
exception of L^. The models also show the magnitude of peak north-south velocity 
y^fx along the principal east-west axis to be approximately equal to the product of 
Lr, the potential vorticity qf^^ in the collar, and a dimensionless number that depends 
weakly on {Dx)\ILr (see Table CI). Since (Av)i is known from the separation of the 
north-south peaks, and V^^y^ can be measured to within ±7 m s^\ the best-fit values 
of Lr = 2400 km and q^^^ = 10.5 x 10^^ s^' were determined simultaneously by 
demanding that the model reproduce the value of as well as the velocity fall-off 
in the outer portion of the collar. 
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Fig. 9. Solid line is far-field po tential vorticity ei°° for model with best-fit parameters given in Table|2] 
Dashed line is q°° determined in lDowling and Ingersoll Jl989l) . 



c. Hollowness of GRS's potential vorticity anomaly 

The models show a uniform potential vorticity anomaly to be inconsistent with 
the north-south velocity along the east- west axis in the GRS 's high-speed collar. In 
particular, an anomaly with uniform potential vorticity cannot simultaneously cap- 
ture the different rates at which the velocity falls-off in the inner and outer portion of 
the collar. However, a model with core potential vorticity ^2 ^■^ f«60% of the collar's 
potential vorticity q'^^^ , is able to capture both fall-off rates to within the uncer- 
tainties. The "hollowness'l^ of the GRS's potential vorticity anomaly is surprising 
because other Jovian vortices such as the White Ovals, as well as other geophysical 

^We define a hollow vortex to be one in which the absolute value of potential vorticity \q\ has a 
minimum in the vortex core. Note that hollowness is not defined in terms of vorticity co; vortices with 
uniform |^| have an |co| minimum in their cores. 
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vortices such as the Antarctic Stratospheric Polar Vortex and Gulf Stream Rings, have 
uniform \q\ or a. \q\ maximum in their cores. Furthermore, hollow vortices are some- 



times unstable (lMarcuslll990l) . Initial- value simulations show that hollow vortices 
re-distribute their q on an advective time scale so that the final \q\ is uniform or has 
a maximum in the core. This raises the question of how a hollow GRS is stabilized 
over time scales longer than an advective time scale. 

d. Non-staircase far-field potential vorticity 

The best-fit q°° has a positive jump AQi of magnitude 1.9 x 10^^ s^^ near the 
eastward-going jet stream and a larger, negative jump AQ2 of magnitude 5.6 x 10^^ sT^ 
near the westward-going jet stream. Due to these opposing jumps, <7°° in this region 
does not monotonically increase with y. This is surprising because numerical and 
laboratory experiments (see §4a) predict a monotonically increasing "staircase" pro- 
file, with jumps only at eastward-going jet streams. The best-fit profile dete r mined 



here agrees qualitatively with a profile determined in iDowling and Ingersolll (119891) 
using an independent method. 

e. Aspect ratio of GRS 's potential vorticity anomaly 

The aspect ratio of the GRS 's potential vorticity anomaly is defined to be (D^) 1 / (^v) 1 , 
where {Dx)i and {Dy)i are the principal north-south and east-west diameters of the 
anomaly respectively. (Recall that the shape of the the GRS's q anomaly, and {Dy)i 
in particular, are obtained as output from the iterative method in appendix B.) The 
aspect ratio of the anomaly depends on the ratio of qf^^ to the shear of the ambient 
flow in which the GRS is embedded; a larger qf^^ to ambient shear ratio irnplies a 



rounder vortex while a smaller ratio implies a more elongated vortex (|Marcuslll99Cl) . 
It should be emphasized that, in general, the ambient shear at the location of the GRS 
is not identical to the shear of the far-field flow v°°. Instead, as shown in Fig.|4l the 
ambient shear is determined by the interaction of the GRS with the far-field flow. In 
particular, the middle column of Fig. |4] shows that v^^^ is large and produces a shear 
with half the magnitude and opposite sign to the shear of v°°. Therefore, the effect 
of v™^ is to greatly reduce the ambient shear at the location of the GRS. For the 
best-fit model, the aspect ratio of the anomaly is 2.18. If the mitigating effect of v^^^ 
on the shear is eliminated by setting AQi = 0, with all other parameters, in particular 
{Dx)\, unchanged from their best-fit values, then the GRS's anomaly shrinks in the 
north-south direction (i.e., {Dy)i decreases) so that its aspect ratio is increased by 
«28%. 

The panel in the left column and bottom row of Fig. |4] explains the functional 
dependence of v^^^ on y and why its shear is adverse to the local shear of v°°. The 
panel shows that the effect of deflecting the jet streams and associated isocontours of 
q°° around the GRS is equivalent to placing nearly semi-circular patches of q north 
and south of the GRS. When the isocontours of q°° that are deflected south of the 
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GRS have latitudinal gradient dq°° /dy > 0, the semi-circular patch of q south of 
the GRS produces anticyclonic shear at the latitude of the GRS. Similarly, if the 
isocontours that are deflected north of the GRS have dq°° /dy > 0, the semi-circular 
patch of q north of the GRS produces cyclonic shear at the latitude of the GRS. Thus 
if the eastward-going and westward-going jet streams of \°°, which are deflected 
respectively south and north of the GRS, both had dq" /dy > 0, then the two semi- 
circular patches of vorticity in Fig. |4] would have opposite sign and form a dipole. 
The dipole would create a large net westward flow at the latitude of the GRS, but 
would create little shear (none, if the patches had equal strength) there. However, for 
the best-fit model, the westward-going jet stream has dq°°/dy < and the eastward- 
going jet stream has dq°° /dy > 0. Both semi-circular patches are anticyclonic and the 
result is a large shear that is adverse to the shear of v°°, as shown in Fig. ID 

7. Conclusions and Future Work 

In this paper we have described a technique for determining quantities of dynam- 
ical significance from the observed velocity fields of a long-lived Jovian vortex and 
its neighboring jet streams. Our approach was to model the flow using the simplest 
governing equation and the fewest unknown parameters that would reproduce the ob- 
served velocity to within its observational uncertainties. For the Voyager 1 data, this is 
a nine-parameter model that is an exact steady solution to the 1.5-layer reduced grav- 
ity QG equations. The nine parameters are the local Rossby deformation radius, the q 
in the GRS's high-speed collar, the q in the GRS's core, the east- west diameter of the 
GRS's q anomaly, the east- west diameter of the GRS's core, the size and steepness 
of two jumps in the far-field q, one located near the latitude of the eastward-going 
jet stream to the south of the GRS, and the other located near the westward-going 
jet-stream to its north. We determined "best-fit" values for the nine parameters by 
identifying several "traits" of the observed GRS velocity field and seeking a model 
that reproduced all those traits. 

Perhaps the most surprising result of our study was that the simple model de- 
scribed above was able to reproduce the entire observed velocity field in Fig. 1 to 
within the uncertainties of 7% (that is, 7 m s^^). The success of the model is due, in 
part, to the fact that the GRS must be well-described by the QG equations, and to the 
fact that the model is an exact steady solution of the governing equations. The suc- 
cess is also due to the fact that the chosen traits are robust and in some sense unique 
(e.g., hollowness) to the physics associated with the GRS. Finally, a part of the suc- 
cess of the model is due to the relatively large uncertainties (7%) of the Voyager 1 
velocities compared to more recent data sets (see below). 

Our most important result was to show that the interaction between the GRS 
and its neighboring jet streams determines the shape of the GRS's q anomaly. By 
explicitly computing the interaction, we showed that the effect of the GRS is to bend 
the jet steams (identified by their jumps in q) so that they pass around the GRS, 
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and the effect of the bending of the jet streams is to reduce the zonal shear at the 
location of the GRS. A smaller zonal shear at the location of the GRS compared 
to the q of the GRS implies a smaller east-west to north-south aspect ratio for the 
GRS's q anomaly. The best-fit model has a positive jump at the eastward-going jet 
stream and a larger, negative jump at the westward-going jet stream. The bending of 
these opposing jumps significantly reduces the zonal shear at the GRS, making the 
aspect ratio of the GRS's q anomaly ^^28% smaller (i.e., rounder) than it would be if 
there were no interaction with the jet streams. It is also interesting to note that due 
to the opposing jumps, the far-field q does not monotonically increase from south 
to north, which is contrary to numerical and laboratory experiments that predict a 
monotonically increasing "staircase" profile. 

The GRS's potential vorticity anomaly was found to be "hollow" with core po- 
tential vorticity Ri60% that of the collar; this is curious because hollow vortices are 
generally unstable. The locations of peak velocity magnitude were found to be accu- 
rate markers of the boundary of the GRS's q anomaly, which implies that the area and 
aspect ratio of the anomaly can be inferred directly from the velocity data. On the 
other hand, clouds associated with the GRS are not an accurate marker of the anomaly 
as they differ from the anomaly area by k,50%. This suggests that cloud aspect ra- 
tios, areas, and morphologies should not be used to determine temporal variability of 
Jovian vortices. 

In devising the model, our philosophy was to include no more complexity than 
was required to match the observed velocity to within its uncert ainties. However, 



lower -uncertainty measurements of the velocity field using CIV (|Asay-Davis et al. 



2006|) on observations from Hubble Space Telescope, Cassini, and Galileo, may re- 



quire that the QG approximation be r elaxed in favor of shallow-water. Also, if ther- 
mal observations (|Read et al.l l2006ah are to be accounted for, governing equations 
that permit 3D baroclinic effects will be required. Modeling different data sets would 
show how the best-fit parameter values evolve with time. 

A companion paper to this one shows that the best-fit model is stable and explores 
the stabilizing effects of the hollow GRS-jet stream interaction. Demonstrating sta- 
bility is important because hollow vortices are usually unstable. Finally, there are 
several questions raised by our best-fit model of the GRS that will need to be an- 
swered. How did a hollow GRS form? Why are there no other hollow Jovian vortices 
(for which the velocity has been measured, cf. the current Red Oval and the three 
White Ovals at 33°S, which existed between the mid-1930's and 1998)? One pos- 
sible answer to the second question is that Jovian vortices apart from the GRS lack 
opposing jumps near their neighboring jet streams and the associated reduction in 
shear due to the v ortex-jet stream in teraction. Indeed, a preliminary best-fit model of 
the White Ovals (|Shetty et al.ll2006h does not show opposing jumps near the neigh- 
boring jet streams. 
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APPENDIX A Method for removing spurious velocities and correcting for the 

curvature of cloud trajectories 

The method invo lves two stages of iteration. We start with the velocity from 



Mitchell et all (119811) in which the trajectories are assumed to be straight lines and 
the velocities are assumed to be located mid-way between tie-point pairs (the initial 
and final coordinates of a cloud feature in a pair of images is defined to be a tie-point 
pair). We then spline the irregularly spaced tie-point velocities onto a uniform grid of 
size 0.05° X 0.05°. The first step of the inner loop of the iteration computes, for each 
tie-point pair, the curved trajectory that a passive Lagrangian particle would follow 
beginning at the initial tie-point location ixi,yj) to its final tie-point location (xF^yp), 
using a 5''' -order Runge-Kutta integration. To carry out the integration, the velocity 
field is spline-interpolated from the grid to the off-grid locations where it is required 
by the integrating scheme. The integration creates a set of closely spaced points, 
{xi,yi), i = 1,2,3,- •• ,N, along the trajectory, where {xi,yi) = {xj,y]). In general, 
this trajectory does not end with {x]\i,yN) equal to {xf,yf) as desired. We therefore 
compute a second trajectory {Xi,Yi) starting from the final tie-point location {xf,yf) 
and integrate backward in time using the interpolated velocity from grid points. A 
third trajectory (x,-,y,-) = [{N — /)(x;, j,-) + (/ — l){Xi,Yi)]/{N — 1) is then computed 
as a lineal" interpolation that, by construction, starts at {xi,yj) and ends at {xf,yf). 
Moreover, because the points along each trajectory are close together, the velocity at 
each point {xi,yi) is well-approximated with the temporal, second-order finite differ- 
ence derivative using the nearest neighbor trajectory points. A new velocity at the 
grid points is created from the spline of the velocities along the curved trajectories of 
all of the tie-point pairs (for each trajectory, we use the velocities at eight approxi- 
mately equally spaced points along the trajectory). We then return to the first step of 
the inner loop. We use the original set of tie-point pairs, but the velocity is now the 
updated velocity on the grid. The inner loop is iterated until it converges (typically, 
three iterations). We then compute the residual of each velocity vector, which is de- 
fined to be the magnitude of the difference between the original, uncorrected tie-point 
velocity and the converged velocity interpolated by splines to that location. Velocity 
vectors with residuals that were six times the root-mean-squared value of all of the 
residuals were considered to be spurious, and their tie-points were removed from the 
original data set. Once the spurious points are removed, the outer loop is complete 
and the entire process is repeated starting with the new (diminished) set of tie-points. 
The outer loop was iterated until no more tie points were removed. The Voyager 1 
tie-point set required three iterations of the outer loop and resulted in the removal of 
220 of the original 1 100 points. The root-mean-squared residual of the iterated ve- 
locity is 7 m s^ \ and we use this value as a measure of the uncertainty in the data. 
For comparison, it should be noted that the residual of the Voyager 1 tie-points with- 
out correcting for curvature and without removing spurious tie-points is f» 9 m s^ \ 
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and the residual for the Hubble Space Te lescope data (frorn CIV ) for the GRS using 
the method described here is Rr: 3 m s^^ (lAsay-Davis et alj|2006h . In the high-speed 



collar, the residuals in the vorticity derived by differentiating the Voyager 1 velocity 
are « 35% of the maximum vorticity. 

APPENDIX B Iterative method for computing steady-state solutions of the 
1.5-layer reduced gravity QG equations 

Here we describe an iterative method for computing steady solutions of Eqs.[T]-[2] 
subject to periodic boundary condition^ in x and y. The method seeks solutions that 
consist of a single anomalous patch embedded in a zonal flow, and that are steady 
when viewed in a frame translating with the patch. Such solutions are of the form 
q{x,y,t) = q{x — Ut,y), where Ux is the constant drift velocity of the vortex. Substi- 
tuting for q in Eq.[T]we obtain: 

(v - f/x) ■ = 0, (Bl) 

which implies that isocontours of q and isocontours of ^if + Uy are coincident. It 
is this property that the iterative method exploits to compute uniformly translating 
solutions. As input, the method requires L^, qf^^, {Dx)i for / = 1,2, •■• ,M, and 
AQ;, 5( for / = 1,2. As output, the method provides q'^^ , U, and the shape of each 
vortex patch. Initial guesses must be supplied for the quantities obtained as output. 
The guesses are then iterated, keeping the input quantities fixed, until the total q is 
a uniformly translating solution. The iterative procedure is described below. The 
domain isx G [—Lx,Lx\, y G [— Ly,Ly]. The origin is at the point of intersection of the 
principal axes. 

1. A new \|/ is computed from the current q by inverting the Helmholtz operator 
in Eq.[2l The current q is the sum of q"', the current q'^'^^, and the current <7™^. 

2. A new drift speed U f or the anomaly is computed. The drift speed of the 



anomaly, as derived in iMarcusI (119931) . is given by J^q'^'^^VxdA/ J^q^'^^dA, 



where A is the current area of the anomaly. 

New isocontours of {\\f+Uy) are computed. The isocontours are streamlines of 
the current velocity in the translating frame. Streamlines that extend from the 
western to the eastern boundary of the domain are referred to as open stream- 
lines. Streamlines that are not open are referred to as closed. 



While periodicity is natural in the east-west direction, it is artificially imposed in the north-south 
direction. This is done by embedding the domain of interest (where the velocities are designed to match 
those of Jupiter) into one with 20% larger latitudinal extent. The flow velocities in the northern and 
southern extremities of the enlarged domain do not match those of Jupiter, but smoothly interpolate the 
velocities from the domain of interest to the periodic boundaries. 
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4. A new q is computed. This is done by setting the value of q along an open 
streamline to the value of <7°° at the point on the western boundary through 
which the streamline passes. In other words, ify = s{x) is the equation of an 
open streamline, then q{x,s{x)) = q°°{s{—Lx)) fox x € [—Lx,Lx]. 

5 . A new q^^^ is computed by computing a new boundary f or patches / = 1 , 2 , • • • , M. 
The new boundary is identified as the closed streamline that passes through 

X = {Dx)i/2. Note that if the current patch is reflection symmetric about the 
N-S axis, the value of {Dx)i is conserved. The potential vorticity of each patch 
qf^^ is held fixed. 



The iterations are repeated until U converges to within a desired tolerance, or 
equivalently, until isocontours of {v/ + Uy) and isocontours of q are coincident. For 
all calculations in this paper, the initial guess for the shape of a patch was an ellipse 
with (Dv),- = 0.5 (Dv),. The final shapes are reflection symmetric about the N-S axis, 
but they are not symmetric about the E-W axis. The initial q'^^ and U were set to 
zero. The grid resolution was 0.05° x 0.05°. The equilibria are not sensitive to the 
domain size provided the domain boundaries are at least three deformation radii away 
from the edge of the outermost patch. We note that it would be interesting to explore 
initial guesses that are not reflection symmetric about the N-S axis, to see if asym- 
metry persists for the fina l solution. Indeed, rece nt low-uncertainty measurements 
of the GRS velocity field (|Asay-Davis et all 120060 show asymmetry about the N-S 
axis. For the Voyager 1 data set however, any asymmetry is much smaller than the 
uncertainties, so asymmetric models are deferred to future work. 



APPENDIX C Sensitivity of model traits to model parameters 



Here we quantify the sensitivity of a model trait to small changes in a model 
parameter. The results justify the methodology used in §5 to determine the best-fit 
parameter values. Consider a trait of the N-S velocity along the E-W axis, say the 
peak magnitude V^^^. From dimensional analysis it is rigorous to write: 

(CI) 

where F is a dimensionless function of seven dimensionless arguments (note that 
Eq. ICll is completely general if the value of V^^-^ is independent of v°° as is suggested 
by decoupling; otherwise, and in particular, for any trait of the E-W velocity along 
the N-S axis, the function F would have to include arguments of the dimensionless 
scalars that parametrize v°°). The sensitivity of V^^^ to changes in a particular pa- 
rameter, say Lr, was determined by computing the value of V^^^ for ^ change in 
of ±5% around its best-fit value with all other parameters fixed at their best-fit val- 
ues, and then using a finite difference scheme to construct the dimensionless partial 
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derivative (L^/V^^^) (3y^^/3L,.) = dlnV^^,,/dlnL,.. Dimensionless partial deriva- 
tives computed for other traits are listed in Table ICll We consider a trait to be in- 
sensitive to any parameter for which the absolute value of its dimensionless partial 
derivative is much less than unity. Note that the results are consistent with Table 1. 

The partial derivatives are not independent. For example, four of the parameters 
in Eq. lCll have dimensions of inverse time (and we write them as x,, / = 1,2,3,4), and 
four have dimensions of length (and we write them as X;, /= 1,2,3,4). Differentiation 
of Eg. ICT] yields the following constraints: 

4 

Y,dlnV^i/dlnXi = l, (C2) 

i=l 

and 

4 

"^dlnV^/dlnXi^l. (C3) 

In general, a trait L that has dimensions of length, such as the width of the N-S 
jet, and which depends only on the parameters in Eq. ICll must satisfy the following 
constraints: 

4 

J^dlnL/dlnZi = 0, (C4) 

and 

4 

Y^dlnL/dlnXi = \. (C5) 

i=i 

Table ICll shows that all traits with the exception of {Dy)i satisfy the constraints. The 
reason {Dy)i does not satisfy the constraints is that it is a trait of the E-W velocity 
and therefore also depends on parameters associated with the far-field flow v~. 

The uncertainties in the best-fit parameter values may be quantified as follows. 
The L2 norm difference between the best-fit velocity and the velocity in Fig. 1 is 
computed. The L2 norm difference is then recomputed with all parameters fixed at 
their best-fit values with the exception of parameter (say). A curve of the L2 norm 
difference as a function of is then computed. By construction, the curve has a 
minimum at the best-fit value of L^. The width of the curve at half-minimum is iden- 
tified as the uncertainty in L^. Since measurements of the GRS velocity using CIV 
have much lower uncert ainties than the Voyager velocity and will soon be available 



(lAsay-Davis et al.ll2006h . we did not deem it useful to compute parameter uncertain- 



ties for the analyses in this paper. 
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Table CI. Sensitivity of model traits to model parameters. Each entry is (parameter/trait) x 
3(trait)/3(parameter) = 3Zn(trait)/9/n(parameter). The partial derivatives were evaluated at the best- 
fit parameter values in Table|2] The fuU-width-at-half-maximum (W) for the northward-going jet along 
the principal east-west axis is broken into two pieces W = E + 1. £ (or /) is the distance along the 
east-west axis between the location where the north-south velocity reaches its maximum value of VjJ^^, 
and the location in the outer (or inner) portion of the high-speed collar where the north-south velocity 
falls to Since the models are reflection symmetric about the principal north-south axis (due to 

the reflection symmetric initial guess, see appendix B), the northward-going and southward-going jets 
have identical values of Vj^^, /, and E. 



Model trait | / Model parameter 




L, 


-7^ 


1^' 




AQi 


AQ2 


8i 


1 


Peak N-S velocity along E-W axis, V^^^ 


0.3 


1.0 


1.1 


0.0 


-0.2 


0.1 


-0.1 


0.0 


c 


Distance between N-S peaks along E-W axis 


1.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


c 


Exterior width of N-S jets at half maximum, E 


0.1 


1.0 


0.1 


0.0 


-0.1 


0.0 


0.0 


0.0 


c 


Interior width of N-S jets at half maximum, / 


1.1 


0.5 


0.0 


-0.1 


-0.7 


0.1 


-0.1 


0.0 


c 


N-S diameter of GRS's potential vorticity, {Dy)\ 


0.6 


1.1 


0.4 


-0.1 


-0.3 


0.1 


-0.3 


0.0 
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